Doped graphane: a prototype high-T c electron-phonon superconductor 
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We show by first-principles calculations that p-doped graphane is a conventional superconductor 
with a critical temperature (T c ) above the boiling point of liquid nitrogen. The unique strength of 
the chemical bonds between carbon atoms and the large density of electronic states at the Fermi 
energy arising from the reduced dimensionality synergetically push T c above 90K, and give rise 
to large Kohn anomalies in the optical phonon dispersions. As evidence of graphane was recently 
reported, and doping of related materials such as graphene, diamond and carbon nanostructures is 
well established, superconducting graphane may be feasible. 



The discovery of superconductors such as magnesium 
diborideQ and iron pnictidesjl, Q opened new hori- 
zons in the landscape of superconductivity research, fu- 
eling renewed interest in the quest for high-temperature 
superconductivity in materials other than the copper 
oxidcs0,Q. The critical temperature, T c , reflects the en- 
ergy scale of the quantum-mechanical interactions driv- 
ing the electron condensation into the superconducting 
stateQ. In high-Tc copper oxidesQ the nature of the 
interaction leading to superconductivity is still under 
debate Q, yet it is generally accepted that Coulomb ex- 
change and correlation effects, with energy scales around 
few hundred meVs, play an important r ole@, E3- In 
contrast, in conventional superconductors the pairing is 
known to be driven by the interaction between electrons 
and lattice vibrations, with an associated energy scale of 
only a few ten meVsjll|. Due to the order-of- magnitude 
difference between such energy scales, it is generally as- 
sumed that conventional superconductors cannot exhibit 
T c as high as copper oxides [H, Here, we report first- 
principles calculations showing that p-doped graphane 
would make a conventional superconductor with T c well 
above the boiling point of liquid nitrogen. 

The Bardeen-Cooper-Schrieffer (BCS) theory^ de- 
fines the basic theoretical framework to understand con- 
ventional superconductivity. I ts g eneralization, known as 
the Migdal-Eliashbcrg theory [12j. incorporating the lat- 
tice dynamics, provides a predictive computational tool. 
Within BCS, T c is given by[ll[: 



N F V 



(1) 



where fee is the Boltzmann constant, Hujq a character- 
istic phonon energy, Np the electronic density of states 
(EDOS) at the Fermi Energy, Ep, V an effective pairing 
potential resulting from the net balance between the at- 
tractive electron-phonon couplin g (E PC) and the repul- 
sive electron-electron interactionjllj. Even though the 
original BCS formula for T c is now replaced by more 
refined expressions such as, e.g., the modified McMil- 
lan equation (Hi]. EqQ] still proves useful for discussing 



trends. EqQ] indicates that one could maximize T c by in- 
creasing the materials parameters ojq, Np, V. However, 
these are str ong ly intertwined, making such optimization 
complex (T3L Il4j . Here, we propose a simple procedure, 
based on EqQ] to design a high-T c superconductor. 

Let us first consider the conventional superconductor 
with the highest T c , MgB 2 (T c = 39K)Q. For simplicity, 
we neglect multi-band and anisotropy effects, which were 
the object of detailed investigations 15 - lj| . In MgB 2 the 
EPC contribution to V is large (~1.4eV, from A = NpV, 
using Np = 0.7states/cell/eV, and A ~ l[3|) because 
the states with energy close to Ep (those which conden- 
sate in the superconducting state ll|) are of a character, 
i.e. derive from bonding combinations of planar B sn 2 hy- 
brids localised around the middle of B-B bonds [15l4l9j . 
These electronic states are significantly affected by the B- 
B bond length variation associated with bond-stretching 



Ei g phonons[16|, |2Cj, resulting into a large EPC contri- 
bution to y. At the same time, the Ei g phonon energy 
is large (~ 60meV[l6j|). due to the small B mass, lead- 
ing to a large ojq in EqQ] Furthermore, MgB 2 is a metal 
with a significant EDOS at ff F (~0.7states/cell/eV[l6|). 
These three factors cooperate in Eq Q] to establish a su- 
perconducting state with T c = 39k |i5I4i9| . However, 
many attempts to improve upon MgB 2 , by investigating 
related materials, only met limited success [2l[, with the 
experimental T c never higher than MgB 2 . 

We thus search for an alternative material having at 
least some of the desirable features of MgB 2 , i.e.(i)tr 
electrons at the Fermi surface, (ii) large bond-stretching 
phonon frequencies, and (iii)large EDOS at Ep. We 
note that the first two requirements are both met by 
B-dopcd diamond, a conventional BCS superconductor 
with T c = 4K[22j], where a small hole- like Fermi surface 
appears around the top of the valence band(23|. The 
electronic states at Ep have a character deriving from 
the bonding combination of tetrahedral C sp 3 hybrids, 
bearing some analogy to MgB 2 . As these a states are lo- 
calized in the middle of the C-C bonds, they couple con- 
siderably to bond-stretching phonons [23j, |24|] , resulting in 
a large EPC contribution to V, even superior to MgB 2 
(~ 3eV, from A = N F V, using Np = O.lstates/cell/eV, 
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FIG. 1: (Color online) (a) ED OS per carbon atom of Id (nan- 
otube; diamond nanowire), 2d (graphene; graphane) and 3d 
(diamond) systems. With the exception of graphene, with 
linear dispersions, the EDOS is proportional to E^ 1 ^ 2 in Id, 
a step-like function in 2d, and E 1 ' 2 in 3d. The step-like 
EDOS in graphane implies that Nf is large even at low dop- 
ing. (b)EDOS of pristine (solid black line) and 12.5% p-doped 
graphane (dashed red line). The top of the valence band 
is set as zero, and Ef =-0-96eV (green line). The EDOS at 
Ef is similar in the two models (0.26 states/eV/cell in rigid- 
band and 0.27 states/eV/cell in supercell).(c)Band structure 
of pristine (solid black line) and 12.5% p-doped graphane 
(dashed red line), (inset) Ball-and-stick 2x2 supercell with 
one substitutional B(top and side views) 



and A ~ 0.3[2^|)(2l In addition, the light C atoms 
have high energy optical phonons (~ 130meV, even after 
softening induced by the large EPC[24[ [25j]). However 
in B-dopcd diamond the EDOS at Ep is rather small 
(-0.1states/ccll/cV for 2% doping^ [H]). This com- 
promises T c . Thus, while B-dopcd diamond shares some 



of the desirable features of MgB2 , its 3-dimensional (3d) 
nature implies that the EDOS in proximity of the valence 
band scales as^ E 1 / 2 (with E measured from the valence 
band edge) [2^1, FigOJa). Then, the number of carriers 
available for the superconducting state remains relatively 
small even for large doping. Superconducting diamond 
is thus a 3d analogue of MgB9[23l.l24[. 

This leads to the question of what would happen in 
a hypothetical B-dopcd diamond structure with reduced 
dimensionality, such as a thin film or a nanowire, where 
the EDOS can be significantly enhanced by quantum con- 
finement. Indeed, the EDOS of a two-dimensional semi- 
conductor goes as~ 9{E) {9 being the step function) (29| . 
hence the number of available carriers can be large, even 
at low doping. In order to estimate the expected EDOS 
increase in a diamond thin film it is helpful to consider 
a simple parabolic band model. For 2% B doping, bulk 
diamond has Np=Q.l states/eV/cell at Ep. A 0.5nm 
thick diamond film with the same doping would have 
Np ^0.5 states/eV/cell. Such an EDOS increase would 
significantly enhance T c . Using the electron- phonon po- 
tential and the phonon frequency of bulk diamond, Eqfl] 
gives that a 0.5nm film would supcrconduct at T c ~ 80K. 
However the question remains whether it is possible to 
synthesize an atomically thin diamond film. 

Recent work on graphene and its derivatives points to a 
positive answer. Soon after the discovery of graphene [3oT| 
several works considered how to functionalise and chem- 
ically modify this novel 2d material [31r!36j . In partic- 
ular, it was proposed that fully hydrogenatcd graphene 
(graphane) would be stable 37] . The main difference be- 
tween graphene and graphane is that, while the former 
is fully sp 2 bonded, the latter is sp 3 , as diamond (37j. 
Recently, some experimental evidence of graphane was 
reported [38j. Since graphane is the 2d counterpart of 
diamond, our scaling arguments immediately point to 
doped graphane as a potential high-T c superconduc- 
tor. Doping could be achieved by gating, including us- 
ing an electrolyte gate, or by charge-transfer, as done 
in graphene 31-34, HiJ H(J. Substitutional doping of 



graphene was also reported, up to ~ 10 cm 2 [351. 141 



We thus perform density functional perturbation the- 
ory (DFPT) calculations of EPC and superconductivity 
in doped graphane within the framework of the Migdal- 
Eliashberg theory |1 211 and the local density approxima- 
tion (LDA)[42|, |43| . By analogy with B doped dia- 
mond, we consider p-doping. This is simulated using 
the rigid-band approximation [44j . Figfljb) shows that 
the calculated EDOS in p-doped graphane close to the 
valence band maximum follows a step-like behavior, as 
expected for a 2d system. At 3% doping the EDOS is 
0.22 states/eV/cell, compared to 0.13 states/eV/cell in 
bulk diamond, with a factor 1.7 enhancement. Figfljc) 
indicates that the dispersions close to Ep are essentially 
identical for a supercell containing B and for a rigid-band 
model of doped graphane. We expect this to hold also 
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FIG. 2: (Color online) ( a) Phonon dispersion of pristine (solid 
black line) and 1% p-doped graphane (dashed blue lines). 
The C-H stretching modes have higher frequencies (2655- 
2711cm -1 ) and are not shown, (b) Optical modes around the 
zone centre, showing the Kohn Anomalies. The horizontal 
(green) arrows indicate the average Fermi surface diameter. 



for lower doping, where the perturbation to the pristine 
dispersions is smaller. The similarity between these two 
models justifies our use of the rigid-band approximation. 
A supercell calculation with the B dopant explicitly in- 
cluded does not show impurity states inside the gap. 

Figf21Ja,b) report the phonon dispersions of pristine 
and p-doped graphane and Fig(3Ja) the corresponding 
phonon density of states (PDOS). Upon doping, the opti- 
cal zone-centre modes with in-plane C-C stretching soften 
as a result of the inception of Kohn Anomalies [45|. The 
two degenerate TO modes, having planar C-C stretching 
and H atoms moving in-phase with the C atoms, down- 
shift from 1185 to 715cm- 1 (147 to 89meV). This is due 
to the large EPC of planar C-C stretching, which signif- 
icantly affects the sp 3 Tike electronic states at the Fermi 
surface. The two degenerate zone-centre modes, hav- 
ing in-plane C-C stretching and H atoms moving out-of- 
phase with respect to the C atoms, downshift from 1348 
to 1257cm- 1 (167 to 156mcV). The LO mode, with out- 
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FIG. 3: (Color online) (a) PDOS of pristine and doped 
graphane. (b) Eliashberg function in p-doped graphane for in- 
creasing doping. The largest contribution comes from the 
optical modes, similar to diamond [24l]. but also the acous- 
tic phonons couple to holes at the Fermi surface, similar to 
SiC [441] . (c) Contributions to the Eliashberg function arising 
from the TO stretching modes (hashed region), (insets) Ball- 
and-stick representations of two TO modes. The arrows indi- 
cate the in-plane C-C stretching motions (carbons are shown 
in grey, hydrogens in white) 



of-plane C-C stretching, does not couple to the electrons 
due to the different parity of potential and wavefunctions, 
resulting into a vanishing EPC. The two degenerate opti- 
cal modes corresponding to the shear motion of the C and 
H planes (at~1133cm- 1 ) and the C-H stretching modes 
(2 modes at 2663 and 271 lcm -1 ) do not undergo soften- 
ing upon doping. This is consistent with the electronic 
states associated with the C-H bonds having little weight 
at Ep, hence a small EPC. 

The softening of modes with a large C-C stretch- 
ing component is similar to that reported in B-doped 
diamond 124 
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In particular, the region of recip- 
rocal space where the phonon softening is observed 
matches the diameter, 2kp, of the hole Fermi surface 
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FIG. 4: (Color online) (a)EPC of graphane as a function of 
doping, calculated using the standard DFPT formalism [HH]: 
the Brillouin zone is sampled with an electron grid up to 
300 x 300 x 1, smearing from 50 to 270 meV, and phonon 
grid of 100 x 100 x 1. For comparison, we plot literature 
values for MgB2 (solid red line[24|]), CaC6 (dashed green 
line[52T|1. and diamond (solid black line[2^|; triangles[25| ) . 
More sophisticated calculations taking explicitly into account 
a substitutional dopant, such as B, could slightly change 
the EPCs[25|, [2S|. However, in B-doped diamond a rigid- 
band model provides a lower EPC and a lower bound to 
T c [23|.(b)T c calculated using the modified McMillan formula 
and a Coulomb pseudopotential /i* = 0.13[S§]. The left- 
side hashed region indicates doping below the estimated MIT, 
where our formalism applies only to charge-transfer or gate- 
induced doping. Above MIT, it applies to substitutional dop- 
ing as well. We use the isotropic Eliashberg formalism [4Sj. A 
more sophisticated description based on the fully anisotropic 
Eliashberg theory is expected to increase T c [l5|, [H^]. For com- 
parison we also show T c of MgB2 (solid red line, T c = 39K[3]), 
CaC6 (dashed green line, T c = 11.5K[H3l), and diamond (solid 
black line, T c = 4K at~3% B[H; UK at~ 7% B[H 



around the T point, this being a typical signature of 
the Kohn effect 45j . The calculated phonon softening of 
the TO C-C stretching modes (~ 58mcV or~ 470cm" 1 ) 
is significantly larger than in other materials, as typi- 
cal Kohn anomalies range from ~ 5meV (graphite and 
graphenefii]) to~ lOmeV (TaC [13). In the case of B- 



doped diamond the phonon softening takes place through 
the creation of a non-dispersive defect branch associated 
with the B dopant [2^|. A similar effect could happen in 
B-doped graphane, but we expect the magnitude of the 
doping-induced softening to be reasonably well described 
within our rigid-band model. Also, more sophisticated 
calculations, taking B explicitl y in to account [25J, [28J or 
with non-adiabatic corrections 48 1, may slightly revise 



the softening. Nevertheless, such a large softening stands 
out as a qualitative effect. 

Figure [3fb) plots the Eliashberg spectral function j49|, 
which shows the relative contribution of different modes 
to the superconducting pairing [49^: 



a 2 F (uj) = i^cj ql ,A ql ,(5 (l 



(2) 



where A qt , is the EPC for a phonon mode v with momen- 
tum q and frequency w qjy , and S is the Dirac delta (we 
used a Gaussian of width 2meV for this purpose). We 
get that the TO in-plane C-C bond-stretching phonons 
with C and H atoms moving in-phase (see FigJ3Jc)) have 
the largest EPC, due to the a character of the elec- 
tronic states at Ep and the large C displacements asso- 
ciated with these modes. This is similar to B-doped bulk 
diamond (23- 25. 28| and validates our hypothesis that p- 
doped graphane can be regarded as an atomically thin 
diamond film, exhibiting similar EPC and vibrational fre- 
quencies, but larger EDOS at Ep. We note that the in- 
plane C-C bond-stretching phonons, with C and H atoms 
moving out-of-phase, do not contribute to the EPC. This 
happens because, upon softening, the four C-C planar 
stretching modes hybridize in such a way that those at 
715cm" 1 carry an increased weight on the C atoms, while 
the opposite happens for the two modes at 1257cm" 1 . 

Figure 0|a) plots the EPC as a function of doping, and 
Fig. IDJb) the corresponding T c . We find that T c exceeds 
the boiling point of liquid nitrogen, and falls within the 
same T c range as copper oxides [501] . Due to the rela- 
tively constant EDOS below the top of the valence band 
[cf. FigfTJb)], T c is rather insensitive to doping. This is 
important for the practical realization of superconduct- 
ing graphane. Our results should be valid throughout 
the entire doping range considered here in the case of 
gate- or charge transfer-induced doping, since in these 
cases the holes are delocalized and doped graphane is 
in the metallic regime. On the other hand, for substi- 
tutional doping we expect our results to be valid only 
beyond the Mott metal-to- insulator transition (MIT). In 
absence of experimental MIT measurements in graphane, 
we estimate the critical doping concentration, n c , using 
the following argument. In 3d the MIT occurs when the 
impurity wavefunctions are close enough that their over- 
lap is significant [5S]. For many materials a^n^ 3 ~ 0.26, 
an being the radius of the ground-state wavefunction of 
an hydrogenic donor 56]. The radius can be calculated as 
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an = e/m*ao/2, a being the Bohr radius, e the dielec- 
tric constant, and m* the effective massjHril]. In diamond 
an ~ 4A and n c ~ 4 • 20cm -3 [56j|, therefore the average 
separation between nearest neighbor B atoms is~ 15A. 
In the case of graphane we use a similar criterion, re- 
placing the 3d hydrogenic impurity with a 2d one. The 
ground-state hy drog enic wavefuction in 2d has a radius 



,2(/ 
'H 



e/m*ao/2[57j. Using the dielectric constant and 
hole effective mass of diamond (e = 5.7; m* = 0.74) we 
find a^f = an/2 ~ 2A. Thus, the average separation be- 
tween nearest neighbor B atoms at the MIT is ~ 7.5A 
and the corresponding doping can be estimated as 5% 
B (1 B every 20 C atoms) or 2 ■ 10 14 holes-cm" 2 . This 
could be feasible, considering that substitutional doping 
in graphene was reported up to 5%[41|. 

The calculated high-T c for p-doped graphane bears 
consequences both for fundamental science and appli- 
cations. One could envision hybrid superconducting- 
semiconducting circuits directly patterned through litho- 
graphic techniques, graphane-based Josephson junctions 
for nanoscale magnetic sensing, and ultimately an ideal 
workbench for exploring the physics of the superconduct- 
ing state in two dimensions [58j. The superconducting 
phase transition in graphane could also be controlled 



by gating [34|, |59(. A high-T c superconductor with gate- 
controllable T c could lead to novel switching mechanisms 
in nanoscale field-effect transistors. Furthermore, the 
discovery of an electron-phonon superconductor with T c 
above liquid nitrogen would mean that(i) there are no 
fundamental reasons to believe that BCS superconduc- 
tors cannot have T c >40K (MgB 2 ), and (ii)high-T c su- 
perconductivity does not take place exclusively in the 
copper oxides. In particular, our calculations indicate 
that at least one material could exist where a very strong 
EPC leads to T c in the copper oxide range without trig- 
gering a lattice instability. The superconducting phase 
transition in systems with reduced dimensionality has 
been the subject of numerous theoretical studies [60j, [61 1. 
Quantum fluctuations could destroy the superconducting 
order in 2d[62|. However, recent experimental evidence 
suggests that this is not necessarily the case 58. 59l 63j|. 
In particular, for thin Pb it was reported that the super- 
conducting state is robust down to two atomic layers 58j . 
Since our proposed mechanism of superconductivity in 
doped graphane is BCS-like, as in Pb, there should be no 
fundamental limits to prevent the realization of high-T c 
superconductivity in graphane. 

It is immediate to extend the present study to diamond 
nanowires, which have been the subject of intense inves- 
tigations in the past few years 64j. For a Id system the 
EDOS near a band edge has a van Hove singularity going 
as~ £J -1 / 2 [29[. We can assume phonon energies and EPC 
to be similar to bulk diamond and graphane. Then EqJT] 
would yield T c as high as~ 150K for a lnm nanowire (see 
EDOS in Figfjja)). The possibility of achieving T c higher 
than copper oxides by exploiting dimensionality deserves 



further investigation. Our work suggests that p-doped 
diamond nanostructures have an intriguing potential for 
high-r c BCS-like superconductivity. 
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